Crossover behavior for long reptating polymers 
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convenient to introduce a small driving external field e, 
applied in a direction tilted by 45° degrees with respect 
of the axes of the lattice: following previous work |$ we 
assign a rate B = exp(e/2) for moves of reptons in the 
direction of the field, while moves in the opposite direc- 
tion occur with a rate B~ l . Here e is a dimensionless unit 
for the strength of the driving field. We focus here on the 
properties in the limiting regime of small e, although the 
DMRG technique is not restricted to this regime. 
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We analyze the Rubinstein-Duke model for polymer reptation by means of density matrix renor- 
malization techniques. We find a crossover behavior for a series of quantities as function of the 
polymer length. The crossover length may become very large if the mobility of end groups is small 
compared to that of the internal reptons. Our results offer an explanation to a controversy between 
theory, experiments and simulations on the leading and subleading scaling behavior of the polymer 
renewal time and diffusion constant. 
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The study of the dynamical properties of polymers is 
a field of great interest, because of important applica- 
tions ranging from material science to biophysics. The 
process of reptation, i.e. the motion of a polymer along 
its own contour by the diffusion of stored length, is gen- 
erally believed to be one of the most important mecha- 
nisms for polymer dynamics pi§|. The simplest model 
for reptation is that introduced by Rubinstein fl and 
later extended by Duke 0] to include the effect of a driv- 
ing field. In spite of its simplicity the Rubinstein-Duke 
(RD) model contains the essential physics of reptation 
and compares well with experiments Q. As very little 
exact results are available for the model one has to rely 
on numerical techniques to investigate its properties. 

In this Letter we present results on the scaling 
behavior of the polymer renewal time and diffusion 
constant, which are obtained by the Density-Matrix 
Renormalization-Group (DMRG) technique ||. DMRG 
allows to compute stationary state properties for rather 
long polymers with unprecedented accuracy. We find a 
crossover behavior with two different regimes: for an in- 
termediate range of lengths various quantities scale as a 
function of the polymer length with effective exponents 
that differ from the asymptotic ones. We also show that 
by tuning an appropriate parameter it is possible to in- 
crease or decrease considerably the crossover length, i.e. 
the characteristic length which separates the two regimes. 
The findings help to clarify controversial results between 
theory, simulations and experiments for the renewal time 
and diffusion constant. Crossover behavior for single 
polymer reptation has been suggested earlier on various 
occasions . Our results are an articulation and proof 
of these suggestions. Moreover they enable us to pin- 
point the crossover region precisely and to estimate the 
effective exponents for arbitrarily long polymers. 

We consider here a d-dimensional version of the RD 
model on a hypercubic lattice (see Fig.[l](a)). The poly- 
mer is divided in N segments, or reptons, of the size of 
the order of the persistence length and each lattice site 
can accommodate an unlimited number of them. It is 




FIG. 1. (a) Configuration of a chain with N = 8 reptons 
embedded in a two-dimensional lattice, (b) One-dimensional 
projection of the configuration of (a) along the direction of the 
applied field identified by the sequence of relative projected 
coordinates y = {1,0, 1,1,— 1,1,0}. The vertical arrows rep- 
resent the allowed moves for the reptons. 

The stationary state properties of the system can be 
derived from the solution of the Master equation 



dP(y,t) 
dt 



[W(y\y')P(y', t) - W(y'\y)P(y, t)\ 



y' 



= Y J H yy .P{y',t) 
v' 



(1) 



in the limit t — > oo. Here P(y,t) indicates the probabil- 
ity of finding the polymer in a configuration y at time 
t and W(y'\y) is a transition rate per unit of time from 
a configuration y to a configuration y' . The matrix H 
contains both the gain and loss terms and is stochastic 
in the sense that the sum over all columns vanishes, as 
required from the conservation of probability. 
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Since the transition probabilities depend only on 
the projected coordinate along the field direction the 
RD model becomes essentially one dimensional (see 
Fig.[l](b)). The relative coordinates between neighboring 
reptons can assume only three values y = {— 1,0, +1}, 
therefore for a chain of N reptons there are 3 W_1 possible 
configurations. One should distinguish between moves 
for internal and end reptons. In terms of the y coor- 
dinates the moves are: (a) Exchange of 0's and l's for 
internal reptons, i.e. ±1, <-> 0, ±1, (b) end repton con- 
tractions ±l—>0 and (c) end repton stretchings — > ±1. 

The only effect of the dimensionality appears on the 
rates for the moves (c), which are dB (dB^ 1 ) for moves 
in the direction of (opposite to) the field, as the end rep- 
ton can move to d unoccupied new sites (here dimension 
stands basically for the lattice coordination number 
Rates for the moves of type (a) and (b) are not affected 
by d. In general the dimensionality is believed to have 
no influence on the dependence of the properties on N, 
as far as the exponents of the asymptotic behavior are 
concerned. On a quantitave level, as we shall see, the in- 
fluence of d is however large. Higher dimensions stretch 
the polymer by giving a larger rate to the stretching pro- 
cesses (c). Not only integer values of d are physically 
relevant, since the dimension enters in the RD model 
as the ratio between end-point stretching and end-point 
contraction. This allows to consider also values d < 1. 
The limit of small d corresponds to the case where the 
motions of type (c) are suppressed. We will show that 
d strongly influences the crossover behavior and the ap- 
pearance of effective exponents. 

In principle a direct numerical diagonalization of the 
matrix H yields, for instance, the polymer drift velocity 
and the relaxation time [pi . In practice such calculations 
are restricted to polymers of rather short lengths (N < 20 
HU) as ^e dimensionality of the configurational space 
grows exponentially with the number of reptons N. This 
problem can be overcome by the DMRG technique || 
that we apply here to the RD model. DMRG has been 
used for quite some time now and it is known for its accu- 
racy and for the possibility of treating large systems using 
an efficient truncated basis set. Although originally in- 
troduced for hermitian matrices jl(| , the DMRG method 
has also been applied to a series of non-hermitian prob- 
lems jll]]. In the latter case the results are known to be 
less accurate than in the hermitian case, however since we 
are primarily interested on the small field limit the ma- 
trices to be diagonalized are only weakly non-hermitian, 
therefore DMRG is expected to perform well in this case. 

Renewal time - We set first e = 0. In this limit the 
matrix H is symmetric since moves in the direction and 
opposite to the field are equally probable. The longest 
relaxation time r, also known as polymer renewal time, 
is a quantity also accessible to experiments since it is 
related to viscosity. Theoretical arguments predict a 
scaling behavior as function of the chain length of the 



type r ~ N z , with z = 3 0. Various experiments 
yielded numerical values of the exponents systematically 
higher z « 3.2 — 3.4 (see, for instance, [0-|3[|^] and ref- 
erences therein). This apparent disagreement has gener- 
ated some discussions in the past years and quite some 
effort has been devoted to reconcile theory with exper- 
iments. Doi argued that the discrepancy is due to 
finite size effects and proposed the following expression 



N 




1 — \ Tr 



(2) 



with No a characteristic length such that only for N ^> 
Nq the right asymptotic behavior can be observed. Ru- 
binstein H showed that a numerical calculation in the 
RD model (for chains up to N as 100 reptons) yielded 
r ~ TV 3,4 , however the asymptotic behavior N 3 was not 
observed. 
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FIG. 2. (a) Plot of z N as function of N~ 1/2 for d = 1, 2, 3; 
the dashed line is a cubic fit in powers of N~ 1/2 for the d = 3 
case and the dotted-dashed line is obtained from Eq. ^. (b) 
As in (a) for d = 0.10, d = 0.25 and d = 0.50. Here and in 
all the other figures error bars, unless explicitely shown, are 
smaller than symbol sizes. 

The renewal time is the inverse gap of the matrix H, 
a quantity which is easily accessible in the DMRG calcu- 
lation. To estimate the exponent z we considered: 



zn 



hiTjv + i - hiTjy-1 
ln(JV + l) - bi(N - 1) 



(3) 



which converges to z in the limit N — ► oo. 

In Fig. ^(a) we show z^ for d = 1, 2, 3 plotted as a 
function of 1/VN since this type of correction-to-scaling 
term is predicted by Eq.||. Using a cubic fit in powers of 
1/y/N we find z = 3.00(2) for d = 3, z = 3.002(4) for 
d = 2 and z = 2.99(1) for d = 1, in agreement with the 
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theoretical arguments leading to a renewal time scaling 
as r ~ TV 3 [Q . Figure ||(a) shows also the fitting curve 
for d — 3 (dashed) and the prediction from Doi's theory 
obtained by substituting Eq.(||) into Eq.(||) (dot-dashed). 
The latter compares very well with our data in the region 
of large N, but fails to reproduce the change of curvature 
in zn, which is characteristic for the DMRG results. 

While in the case d = 2, 3 the asymptotic value is ap- 
proached monotonically from above, for d = 1 zjv shows 
a non-monotonic behavior. In a plot of In r^r vs. In TV 
the maximum of Zn corresponds to an inflection point, 
where curvature is absent. Therefore in a range of lengths 
around the maximum the numerical data would be fit- 
ted rather well by an effective exponent z c q = maxjy Zjf. 
This is particularly true in the case of data affected by 
small, but non-negligible error bars, as it happens in ex- 
periments or in computer simulations. A monotonic be- 
havior is less harmful since in a lnTjv vs. In AT plot one 
should be able to distinguish always a non- vanishing cur- 
vature, although it may happen that the data are well 
fitted with effective powers in a range of values of N also 
in that case. We refer to the non-monotonic behavior as 
to a crossover effect. To investigate it more into details 
we considered also d < 1 (see Fig. ||(b)). For d = 0.5 
and d = 0.25 again the data show a non-monotonic be- 
havior in A and the turning point shifts to larger A, 
i.e. crossover effects are more pronounced. For d = 0.1, 
the smallest value considered here, and for the range of 
A investigated, z^ increases monotonically. Assuming 
that the asymptotic exponents do not depend on d, the 
turning point will be reached with much longer polymers 
than those considered here. It is also important to point 
out that the values for the effective exponent that we 
estimate from our data z D g rs 3.2 — 3.4 are exactly in 
the range of values which are found experimentally in 
measures of viscosity . Therefore our results strongly 
support the idea that the differences between theory and 
experiments are due to a crossover effect, as pointed out 
by some authors in the past . The advantage of the 
DMRG calculations is here that they provide clear in- 
formation of both the asymptotic and non-asymptotic 
regions. Moreover, in the case of d small, we find a range 
of lengths for which the renewal time scales with a local 
exponent smaller than 3 (see Fig. ||(b)). The crossover 
to the asymptotic regime where z = 3 is approached from 
above (as in Eq. ||) occurs only for extremely long poly- 
mers. We are not aware of any experiments showing this 
type of behavior; possible candidates could be concen- 
trated polymer solutions with bulky endgroups in such a 
way that endpoint strechings are suppressed. 

Diffusion constant - Crossover behavior appears also in 
other quantities, as for instance in the scaling of the dif- 
fusion constant D(N) as function of the polymer length 
A. We calculated D(N) by applying a small field e and 
using the Nernst-Einstein relation p2| 
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For the scaling behavior of D(N) one expects 
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D(N) = 



(4) 



(5) 



where the leading term N~ 2 is by now well-understood 
p[-p|, p^[- 16| and it is considered to be experimentally ver- 
ified |17|. In the RD model also the prefactor happens 
to be known exactly |l^Jl^|: A = 2d + 1. The next to 
leading order term has been investigated as well. By re- 
lating the diffusion constant to the renewal time it was 
predicted |lq| that the correction term would be anoma- 
lous, i.e. 7 = 1/2. This prediction is also supported by 
other theoretical arguments |15|. On the other side ac- 
curate Monte Carlo simulation results, done for d = 1, 
could be best fitted with a power 7 2/3 both for the 
RD model |13( ] and also for another model of polymer 
reptation ]16|. This issue is still unresolved. An expo- 
nent 2 /3 is somewhat surprising since, as also seen above 
for the gap one naturally expects A^ 1 / 2 corrections. 
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FIG. 3. Log-log plot of DN 2 - 1/3 vs. N for d = 1. The 
DMRG data are in good agreement with the exact ones and 
with the Monte Carlo results. Inset: blow up of the region 
for large N where the deviation of the DMRG data from the 
slope 2/3 starts being noticeable. 

In Fig. | we show a plot of g N = D(N)N 2 - l/A 
as function of A on a log-log scale for d = 1. The fig- 
ure shows results from exact diagonalization for small 
lattices, Monte Carlo simulations and DMRG results 
which we extended up to A = 57. DMRG data are in 
good agreement with those obtained from other meth- 
ods and follow rather nicely a slope —2/3 in the plot. 
Only a very close inspection of the region of large- A sys- 
tems reveals (see inset) that this slope is not the correct 
asymptotic one. Similar to the gap, this can be best seen 
from the discrete derivative 
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lN hx(N + 1) - ln(iV - 1) ' ( ' 

A plot of 7jv vs. 1/s/N is shown in Fig. |J for d = 1, 2 
and 3. As for the gap, we note a non-monotonic behavior 
for d = 1, where the maximum of 77V is found at about 
0.68, i.e. very close to the power 7 = 2/3. However 
for sufficiently large TV, 7^ clearly deviates from 2/3 to 
smaller values. For d = 2 and 3, 77V is monotonic in N. 

Like for the renewal time exponent zjf, we fit 7 at with 
a cubic curve containing powers of 1/yN 

3 

i=l 

where 7 and aj are fitting parameters. Extrapolations 
yield 7 = 0.51(1) (d = 3), 7 = 0.51(1) (d = 2) and 
7 = 0.48(3) (rf = 1). These results strongly support a 
correction term with 7 = 1/2. 
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FIG. 4. Plot of the effective exponent 7jv as function of 
l/y/N for d = 1, 2 and 3. The dashed line is the value of the 
correction exponent conjectured on the basis of Monte Carlo 
simulations for d = 1. The dotted line is a guide for the eye. 

In conclusion, the above results show that DRMG is 
powerful technique to investigate the properties of a rep- 
tating polymer. In this Letter we have restricted our- 
selves to two aspects: the renewal time and the diffusion 
coefficient, since these quantities have been mostly de- 
bated in the past. The DMRG calculations also yield a 
host of detailed information about the structure of the 
reptating polymer. In a forthcoming publication we will 
present data on the drift velocity for larger values of e 
and on structural properties of the chain such as the av- 
erages (yi} and (yiVi+x), where i labels the position along 
the chain. Here we have shown that the large finite size 
corrections, characteristic for the reptation process, man- 
ifest themselves as effective exponents for the asymptotic 



behavior of the renewal time and the diffusion coefficient. 
The accurate values for a large set of lengths N make it 
possible to determine the corrections to scaling with great 
precision. In this way we find a natural reconciliation of 
the fairly strong theoretical arguments for a renewal time 
exponent z = 3 and the equally pertinent experimental 
findings of values around z w 3.2 — 3.4. Our analysis also 
shows that log-log plots to determine exponents are haz- 
ardous when such large corrections to scaling are present. 
In particular when the effective exponent shows a station- 
ary point as function of TV, it shows up as an inflection 
point which can be easily mistaken for the asymptotic 
region in case of insufficient data. On the other hand the 
method of differential exponents is a very accurate indi- 
cator of the asymptotic behavior. However differential 
exponents require precision data for many values of the 
length N, which neither in Monte Carlo simulations nor 
in experimental measurements are feasible. 

We are grateful to G.T. Barkema for helpful discus- 
sions. 
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